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Scanning tunneling spectroscopy at very low temperature on homogeneously disordered supercon- 
ducting Titanium Nitride thin films reveals strong spatial inhomogeneities of the superconducting 
gap A in the density of states. Upon increasing disorder, we observe suppression of the super- 
conducting critical temperature Tc towards zero, enhancement of spatial fluctuations in A, and 
growth of the A/Tc ratio. These findings suggest that local superconductivity survives across the 
disorder-driven superconductor-insulator transition. 

PACS numbers: 74.50.-|-r, 74.78.Db, 74.81.-g 



A pioneering idea that in the critical region of the 
superconductor-insulator transition (SIT) the disorder- 
induced inhomogeneous spatial structure of isolated su- 
perconducting droplets develops P, 0] , grew into a new 
paradigm |3i]. Extensive experimental research of criti- 
cally disordered superconducting films revealed a wealth 
of unusual and striking phenomena, including nonmono- 
tonic temperature and magnetic field dependence of the 
resistance 0, 13, S 0|i activated behavior of resistivity 
in the insulating state [J d, @, 0, S| > nonmonotonic mag- 
netic field dependence of the activation temperature, and 
the voltage threshold behavior d, 0, • These features 
find a theoretical explanation based on the concept of 
disorder-induced spatial inhomogeneity in the supercon- 
ducting order parameter 

EETEIil,!!!- Numerical 



simulations confirmed that indeed in the high-disorder 
regime, the homogeneously disordered superconducting 
film breaks up into superconducting islands separated by 
an insulating sea fis', . At the same time the direct 
observation of superconducting islands near the SIT jus- 
tifying the fundamental but yet hypothetical concept of 
the disorder-induced granularity on the firm experimen- 
tal foundation was still lacking. 

In this Letter, we report on the combined low temper- 
ature Scanning Tunneling Spectroscopy (STS) and trans- 
port measurements performed on thin Titanium Nitride 
films on approach to the SIT. The local tunneling density 
of states (LDOS) measured at 50 mK reveals disorder- 
induced spatial fluctuations of the superconducting gap, 
A, with both, standard deviation a to the average gap 
and the gap to the critical temperature ratios, a/ A and 
A/Tc, respectively, increasing towards the transition. 

Our samples were thin TiN films synthesized by atomic 
layer chemical vapor deposition onto a Si/Si02 substrate. 
TiNl was a 3.6 nm thick film deposited at 400° C while 
TiN2 and TiN3 were 5.0 nm thick films deposited at 
350°C. TiN3 was then slightly plasma etched in order 
to reduce its thickness. Electron transmission images 



revealed that the films comprise of the densely-packed 
crystallites with a typical size of 4 to 6nm. The sam- 
ples were patterned into the Hall bridges using conven- 
tional UV lithography and plasma etching. It is worth 
noticing that identically fabricated TiN films undergo the 
disorder- and magnetic field-driven SIT 0, Q. Trans- 
port measurements and STS were carried out during the 
same run in a STM attached to a dilution refrigerator. 
The STM Pt/Ir tip was aligned above one of the free 
500 X 500 /im^ contact pads of the Hall bridge. In order 
to probe the LDOS, the differential conductance of the 
tunnel junction, G(V), was measured by a lock-in ampli- 
fier technique with an alternative voltage of 10 /iV added 
to the direct bias voltage. Resistances were measured by 
acquiring both voltage and current with a low frequency 
lock-in amplifier technique in a four wires configuration. 

Shown in Fig. la are the temperature dependences of 
sheet resistances, R\j, for the three samples, which room 
temperature sheet resistances are listed into Table I. All 
three curves show nonmonotonic behavior with a pro- 
nounced maximum preceding the superconducting tran- 
sition. Such a behavior can be described within the 
framework of the theory of quantum corrections for quasi- 
two-dimensional disordered superconductors. Making 
use of theoretical expressions for corrections 17 , S Esl > 
we can find the superconducting temperature Tc 2u\ (see 
Table I). The calculated fits are presented by solid lines 
in Fig. la. The perfect description of the experimental 
curves achieved by means of formulas derived for ho- 
mogeneously disordered systems, shows that relying on 
the macroscopic parameter R(T), we could have expected 
that our homogeneously disordered films possess a spa- 
tially uniform superconducting state. However, using the 
data obtained by virtue of the local STM probe, we reveal 
that due to mesoscopic fluctuations of disorder, supercon- 
ducting properties show inhomogeneities on mesoscopic 
spatial scales. Typical local tunneling data are shown in 
Fig. lb. In contrast with earlier macroscopic tunneling 
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FIG. 1: (a) Sheet resistance Rfj versus temperature for three 
samples. The solid lines are fits according to localization- 
interaction and superconducting fluctuations corrections. 
The legend of panel (a) describes the two panels, (b) Nor- 
malized differential tunneling conductance measured at T = 
50 mK (dots). Spectra are shifted for clarity. The BCS 
fits (solid lines) were calculated with the following param- 
eters: TiNl - A = 260 /xeV and an effective temperature 
Te// = 0.25 K; TiN2 - A = 225 peV, T^ff = 0.32 K; TiN3 - 
A = 154 fieV, Te// = 0.35 K. 



■ I - ■ *■ 



at 



■1 ■ 



A[nV] 
-280 

■275 

■270 

■265 

-260 

■255 

■250 

■245 

-240 



100 

X [nm] 



G(V), normalized 




1000 1500 2000 

Displacement [nm] 

FIG. 2: Top: The colour map of spatial fluctuations of A 
on TiNl. Inhomogeneities of the superconducting properties 
show up on a scale of a few tens of nanometers. Bottom: 
Spectra measured along a straight line on TiN2. The BCS- 
like LDOS fluctuates symmetrically around the Fermi level. 



studies performed with lithographied junctions on Bis- 
muth films [2l[ , our local STM measurements systemati- 
cally showed fully gapped shapes for any degree of disor- 
der indicating the absence of quasiparticles at the Fermi 
level. These BCS-like spectra at low energy allowed us to 
extract A values which were significantly reduced as com- 
pared to A''""= = 730 AieV in bulk TIN [23|. Yet, spectra 
taken at different positions on the surface give different 
values of A. This evidences that the superconducting 
state is spatially inhomogeneous. 

TABLE I: Sample characteristics: -R300 - resistance per square 
at 300 K. Tc - critical temperature determined from the quan- 
tum correction fits. A - average value of the superconducting 
gap. a - standard deviation of the superconducting gap. 





^300 




A 


a 


a/A 






kn 


K 


fieV 








TiNl 


2.45 


1.3 


265 


11 


0.04 


2.37 


TiN2 


2.7 


1.0 


220 


13 


0.06 


2.55 


TiN3 


3.5 


0.45 


160 






4.13 



The top panel of Fig. 2 presents the spatial map of su- 



perconducting gap A measured over the 200 x 150 nm^ 
area with a pitch of 3.3 nm. The characteristic scale of 
the inhomogeneity is estimated as a few tens of nanome- 
ters. The measured 2700 spectral gap amplitudes give a 
gaussian distribution with an average value A = 265 fieV 
and a standard deviation cr = 11 fieY. The spatial fluc- 
tuations are straightforwardly seen in the bottom panel 
of Fig. 2 displaying a color map of 256 LDOS spectra 
measured on TiN2 every 11.7 nm along a straight line. 
Note, that spectra are symmetric with respect to voltage 
direction. We find A = 220 ^eV and cr = 13 /.teV. In 
TiN3, probably because of the plasma etching, the sur- 
face did not suit for spatially regular STS measurements. 
Therefore we performed point-contact spectroscopy, the 
STM tip gently touching the surface sample. From the 
30 BCS-like spectra we obtained the magnitudes of the 
gap scattered in the interval from 125 to 215 /ieV, with 
an average A — 160 /leV. For the three samples, every 
spectra measured in the scanning window of 3 x 3 /im^ 
displayed gaps consistent with the evaluated distribution. 
All the data are summarized in Table I. They show an 
unusually large (as compared to the BCS-predicted 1.76 
value) A/ksTc ratio and an increasing relative standard 



3 



deviation cr/A with disorder. 

When deaUng with STM and spatially resolved super- 
conducting inhomogeneities, the question of the tunnel- 
ing junction quality cannot be eluded. However, the lack 
of states at the Fermi level in our measurements is a 
strong indication of good tunneling conditions. For in- 
stance, for tunneling from a damaged metallic surface 
weakly coupled to the film with a proximity-induced su- 
perconducting LDOS, one would have rather observed 
small A/ZcbTc ratios and a significant amount of states at 
the Fermi level. These fully gapped spectra also rule out 
any spurious inelastic scattering of the tunneling carriers. 
Moreover the diffusive length probed by them in the mV 
range is about the thickness of the film itself which makes 
questionable any distinction between the physical prop- 
erties of the surface and of the film. Finally, we did not 
observe any dependence of these spectra upon the tip- 
sample distance which varied only slightly for TiNl and 
TiN2 when the tunneling current setpoint was changed 
over orders of magnitude. All this makes us confident 
that the discussed superconducting inhomogeneities are 
of the intrinsic origin. 

The appearance of spatial inhomogeneities of the or- 
der parameter was indeed predicted by theoretical cal- 
culations [U, and by numerical simulations 15 1 as a 
result of interplay between superconductivity and Ander- 
son localization in a strongly disordered superconductor. 
The recent work [2^ showed that even in the weakly 
disordered systems far from the localization threshold, 
the mesoscopic fluctuations of the conductance near the 
point of the Coulomb suppression of superconductiv- 
ity [iS] also give rise to spatial inhomogeneities of the 
order parameter. We like to emphasize that in all the 
above works it was the short range, i.e. atomic scale, dis- 
order that caused variations in the physical characteris- 
tics of the system on the mesoscopic scale and thus drove 
it into an inhomogeneous superconducting state. Numer- 
ical simulations [15| demonstrated further that increas- 
ing disorder washes out the BCS singularities by pushing 
states at higher energies. We indeed observe smearing 
out of these coherence peaks. However, we can well de- 
scribe this effect with BCS formulas with the effective 
temperature T^f / which is higher than the expected 250 
mK set-up resolution Note at this point that BCS 
expressions do not describe the observed upturn in the 
LDOS above the gap (see Fig. lb) which is likel y d ue to 
Coulomb interactions [17| not considered in ref. [15[ • We 
conjecture that Tg// is a consequence of the spatial mod- 
ulation of the superconducting gap. The regions with 
the reduced A act as traps for normal excitations. Tun- 
neling of a quasiparticle between the regions with the 
reduced gap gives rise to a broadening, F, of the quasi- 
particle energy levels. This broadening is equivalent to 
some effective temperature, T^ff ~ T/kB, and can be es- 
timated as r ~ {2m* L^), where L is the characteris- 
tic spatial scale of the modulation of the superconducting 




FIG. 3: Critical temperature Tc and mean spectral gap A at 
T = 50 mK divided by their bulk value versus the room tem- 
perature sheet resistance. Black bars represent the measured 
gap ranges for each sample (~ 4(t for TiNl and TiN2). Blue 
solid line is a fit according to Finkelstein's model [24|. Red 
solid line is a guide to the eye. Black triangles are the room 
temperature sheet resistances for the "last" superconducting 
sample (SI - -R300 = 4.48 kfi) and "first" insulating sample (II 
- i?300 = 4.60 kf7) from the Ref. [|. 



gap, / ~ exp(— L/^) is the tunneling integral, and m* is 
the quasiparticle mass. Taking for a crude evaluation 
L ~ ^ « 10 nm (the latter is derived from the experimen- 
tal value of the U£per critical field Bc2{0) = 2.8 T [2(i]) 
and m* ~ 2me [27l |. we arrive at T^/f « IK, which is 
amazingly close to experimental findings given very ap- 
proximate character of our estimate. Note that the effect 
of smearing of the BCS peaks in the average density of 
states of inhomogeneous superconductors was discussed 
in the classical paper by Larkin and Ovchinnikov (2^ . 

Now we are in a position to discuss the evolution of 
the superconducting properties of the TiN films with 
increasing disorder, measured by i?3oo, the room tem- 
perature sheet resistance. Fig. 3 presents A/A''""^ and 
Tc/T^""^ versus i?3oo for the samples under this study. 
^buik ^ 0.73 meV and T^''""= = 4.7 K were measured on a 
much less disordered thick film with i?3oo = 27J1 [i^]. In 
addition, the data on the sample SI, with i?3oo — 4.48 kO, 
and on the sample II, with i?3oo — 4.60 kil from the 
Ref. [1] are shown. 

Important inferences can be deduced from this figure: 
(i) Increasing disorder suppresses Tc towards its com- 
plete vanish. This behavior is well fitted within 
Finkel'stein's model which describes disorder-enhanced 
Coulomb interaction in homogeneously disordered thin 
films 



2j|. In this model, the reduction of Tc follows 



the relation : In 



= 7 ■ 



^In 



l/-f+r/i-Jr/2 



/2r — yi/^+r/4+y/r/2^ 

where r = i?ne^/(27r^?i) and 7 is the only fitting param- 
eter. We found 7 ~ 6.2 consistent with the previously 
reported data on TiN (7 ~ 6.8) ^S)]. Note that the sam- 
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pies SI and II from Ref. [8| are identically fabricated 
TiN films closest to the SIT from the superconducting 
and insulating sides respectively. 

(ii) The variance in the gap magnitude increases with 
disorder. The enhancement of spatial fluctuations leads 
to a strongly inhomogeneous superconducting state when 
approaching the SIT. This observation supports numeri- 
cal calculations 15| . It further agrees with the conclusion 



of [23| that the mesoscopic fluctuations of the Cooper 
attraction constant in Finkel'stein's theory yield strong 
spatial fluctuations of the local order parameter close to 
the critical disorder. 

(iii) The suppression of Tc advances that of A, yield- 
ing anomalously large A/kBTc ratios (see Table 1) as 
compared to the bulk TiN value l^buik ij^^rj.buik ^ g 
(2^ . This can indicate, should the tendency displayed 
on Fig. 3 persisted till the full suppression of Tc, that 
despite vanishing of the phase stiffness at the critical re- 
sistance, the average superconducting gap remains finite. 
Extrapolating the red line over the insulating side of SIT 
(Fig. 3), one then could have concluded that even in the 
insulating phase where Tc = 0, the average supercon- 
ducting gap still survives. Thus the observed trend of 
the faster decay in Tc than in A, together with the in- 
homogeneous state, offers a strong support to hypothesis 
that the so-called "homogeneously disordered " supercon- 
ducting films near SIT [iiailii can in fact be viewed 

as granular-like superconducting structure or Josephson 
phase suggested by Imry et al |3|]. This, in its turn, im- 
plies that the bosonic mechanism of the superconductor- 
to-insulator transition [29^ becomes relevant. However, 
an alternative theory building a superconducting state 
based on fractal wavefunctions close to the Anderson 
transition has been recently proposed 13|. Consistently 



with our data, the authors predict likewise an inhomoge- 
neous superconductor with strong A/fc^Tc ratios and a 
gapped insulator with no coherence peaks due to pairing 
of electrons on localized states. 

In conclusion, our results demonstrate that the be- 
havior of global characteristics, such as R\j{T) and the 
suppression of by increasing disorder in the TiN 
films near the superconductor-insulator transition are de- 
scribed perfectly well by theories developed for homoge- 
neously disordered films. At the same time the mea- 
surements of the local tunneling density of states clearly 
indicate the existence of the inhomogeneous supercon- 
ducting state. Moreover, the data show that LDOS is 
likely to remain gapped beyond the SIT in the insulating 
phase. This suggests the presence of the regions with the 
finite superconducting gap in the insulating side of this 
transition. Several issues related to our results on local 
tunneling measurements call for further research. First, 
the non-constant LDOS above the gap seen in the lower 
panel of Fig. 1 requires thorough analysis taking into ac- 
count Coulomb interactions. Another important ques- 
tion is the role of spatial modulation of the gap. While 



its effect on the global density of states averaged over 
the whole system was discussed in [2^, the consistent 
theoretical description of the local DOS is still lacking. 
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